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ABSTRACT 

A cylindrical Josephson junction with a spatially dependent Josephson coupling which 
averages to zero is studied in order to model the physics of a disk of d-wave super- 
conductor embedded in a superconducting film of a different symmetry. It is found 
that the system always introduces Josepshon vortices in order to gain energy at the 
junction. The critical current is calculated. It is argued that a recent experiment 
claimed to provide evidence for s-wave superconductivity in YBa2Cus07 may also be 
consistent with d-wave superconductivity. 
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In this note I analyze a theoretical model of an experimental system recently 
proposed as a test for the presence of d-wave superconductivity in the cuprate high- 
T c superconductors [1]. The proposed experimental system is shown schematically 
in figure 1. [The inner region is drawn as a circle, the region used in ref. 1 was 
a hexagon.] The inner region, labelled A, is a c-axis-oriented film of high-T c CuC>2 
superconductor with the Cu-0 bonds running vertically and horizontally as shown. 
In the devices actually made [1] the outer region, B, is another c-axis film of high-T c 
superconductor, with the Cu-0 bonds rotated by 45° with respect to region A, as 
shown. I will also consider the case in which region B is an s-wave superconductor. 
In either case the Josephson coupling alternates in sign as one proceeds around the 
perimeter of region A, and averages to zero. The experiment is to measure the 
critical current of the A-B interface. A nonzero critical current was argued [1] to be 
evidence against d-wave pairing because if one chooses a constant phase a in region 
A and a constant phase 0b in region B, then the Josephson energy of the interface 
vanishes and so there can be no supercurrent across it. This argument is incorrect; 
the difficulty is the assumption of constant, spatially independent phases 0a and 0b- 
I shall show that for the two dimensional geometry of interest, the system introduces 
Josephson vortices to gain energy at the boundary at the expense of introducing 
bulk supercurrents. The physics is similar to that of inductive effects in a SQUID, 
where one balances the Josephson energy against the magnetic field energy. The 
dimensionality is crucial to the argument; the results do not apply to the essentially 
one dimensional "corner junction" SQUID experiments recently argued [2] to provide 
evidence of d-wave superconductivity in YBa 2 Cu 3 07_,5. 

In the remainder of this communication I analyze the interplay between bulk cur- 
rents and the Josephson energy more quantitatively, by minimizing the appropriate 
free energy. The free energy has three contributions: the Josephson energy, the su- 
percurrent energy, and the energy in the magnetic field generated by the supercurrent 
flow. I first study the simplest case, that of a very thin film for which the penetration 
depth may be taken to be infinite and all magnetic effects are negligible. Then I 
consider the general case, and finally discuss the relation to experiments. 

Assume first that the only contributions to the energy are the bulk supercurrents 
(i.e. the kinetic inductance) and the Josephson energy of the A-B interface. The 
energy of a bulk supercurrent involves a superfluid stiffness p which I relate to a bulk 
critical current J c by the argument that at the critical current the phase changes by 
an amount of order 1 in a distance of order the coherence length, £ > i- e - Jc£o = P- 
The energy, E, is then 

E = j o rdrde-JXM A ) 2 + j a rdrde-jz^iv^) 2 
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+Ija ( * d9g M {6) cos[0> = a,0) - <j) B (r = a, 6)] (1) 
Jo 

Here A,B are the superconducting phases in regions A and B respectively, Ij 
is the maximum Josephson current of an infinitesimal region of the interface, and 
gM(#) describes the angular dependence of the Josephson coupling. In the case where 
A is d-wave and B is s-wave, we must on general grounds have J 27r d#gM(#) = 
and gM(0 + 7r /2) = — givi(^)- I call this case M=2. In the case considered in ref. 
(1) we must have J 27r d0g M (0) = and g M (0 + vr/4) = -g M (0); I call this M=4. For 
explicit calculations I take the simplest forms consistent with these considerations [3, 4], 
namely g2(0) = sin 26 or g4(0) = sin 40 . I indicate the modifications arising for the 
alternative case of a piecewise constant gyi(6)- 

The theoretical problem is to minimize the energy in (1) subject to the constraints 
of local current conservation at each point of the A-B interface and fixed total radial 
current. Local current conservation implies (h is the outward normal to the A-B 
interface at angle 6) 

jXn • V<P A = J* ton ■ V0 B (2a) 
= I J g(6) S m[ ( t ) A (a,d)- ( t> B (a,d)] (2b) 

In terms of the mean radial current density at the A-B interface, Io, the constraint 
on the total radial current is 

It r 2w 

Io = j2^Jo d69{6) Sin[0 ( °' 9) ~ ° )] (3) 

To perform the minimization, note that Equation 1 implies that A,B satisfy y 2 — 0. 
Assume that the current 27rl a is introduced symmetrically in the center of region A. 
The problem then has circular symmetry; the general solution to the Laplace equation 
may be written down, the solutions in regions A and B may be related using 2a, and 
the result substituted into (1). 
The result is 

/:; X E + ( 4 ) 



J do 2 
-I r T-9m{&) cos[0(0)] 

where 

(f)(6) = <f) A (a,6) -<f) B (a,6) = 0° + ]T [c n sin n6 + d n cos nO] (5) 

n 

and J = J A J B /(J A + J B ) and I = vrlja/J^o. 

Minimizing (5) with respect to the c n and d n at fixed 0o ensures that all of the 
n > Fourier components of the current conservation equation 2b are satisfied. The 



3 



mean current Io (eq. 3) is then fixed by <f>o, which is the angular average of the phase 
drop across the A-B interface. 

From (4) it is clear that the important parameter is I; for I <g; 1 the Joseph- 
son energy is weak while for I ^> 1 it is strong. Further, I increases linearly with 
the perimeter of the disk, because in two dimensions the energy of bulk superflow 
/ d 2 r(v</>) 2 is scale invariant. Consider first I< 1. Clearly the free energy may be 
expanded in powers of the c n and d n . Expanding to leading nontrivial order and 
minimizing gives 

ij/ sin 20 o 

Io = —^T- (6a) 
and _ 2 

E = ~lM Sin2(p0 (6b) 
Observe that E is minimized at a phase difference of ir/2. The bulk supercurrents at 

0o = are shown in fig. 2, for M = 2. One sees that the solution corresponds to 

2M vortices, with alternating positive and negative circulation, centered at the points 

where gu(d) vanishes. Because I « 1, the "core" of the vortex is much larger than 

the vortex spacing, a. If a mean radial current is imposed, O ^ vr/2 and the centers 

of the vortices shift along the perimeter of the circle in such a way that the sum of 

the imposed radial current and the circulating current about each vortex vanishes at 

the points where gu{9) = 0. The motion of the centers of the vortices is shown in the 

inset to fig. 2. 

In the case of a piecewise constant gu.{Q), i n which g M (# = 1 (0 < # < vr/M), 
gu{0) = —1 (n/M < 9 < 2n/M) etc, one finds that the energy is still proportional 
to sin 2 0o but with a slightly different coefficient, and the qualitative structure of the 
vortices is not changed. 

Note that the energy, eq. 6b, is periodic in O with period ir, unlike a conventional 
junction which would have a 2ti periodicity. A 7r-periodicity is unambiguous evidence 
of unconventional superconductivity and may in principal be measured in a SQUID 
experiment. Now the quantity 0o is the angular average of the phase drop across the 
A-B interface, but from the solution of the Laplace equation that led to eq. (4) I also 
find that at a distance d < a from the center of the disk, the phase has a negligible 
angular dependence and is 

0(d) = 0o - - In a/d (7) 

7T 

Thus for I <sC 7r/ln(a/d) ~ 1 the phase in the center of the disk is essentially O and 
so an experiment in which this phase is controlled would measure a 7r periodicity. 

Now consider the opposite limit, I ^> 1. Specialize first to the case O = n/2; then 
it follows from the minimization equations that all the d n = 0. The Josephson term 
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in eq. 4 is then minimized by a piecewise constant (f)(6), equal to ±7r/2 according 
as g(6) is positive or negative. Now a piecewise constant (f)(6) implies c n ~ 1/n and 
therefore a logarithmically divergent J2n nc n- This divergence comes from the familiar 
1/r behavior of supercurrents outside the core of a vortex. To cut off the logarithm 
one must convert the steps in (f)(6) to smooth crossovers extending over a scale e 
which gives the core size of the vortex. For gu(6) = sin(M#) the cost in Josephson 
energy is ~ Ie 2 , while the superfluid energy is ~ lnl/e. Optimizing gives e ~ I -1 / 2 
and E = — I + Olnl. Thus at (f>o = n/2 and I ^> 1 the system gains essentially the 
maximum Josephson energy and the vortices depicted in fig. 2 are pinched at the 
A-B interface into angular regions of size I-V 2 . For a piecewise constant g, the core 
size would be e ~ I -1 and so the coefficient of the In would change by a factor of two. 

To treat general values of 0o a t I 3> 1 one must write (f)(6) as a sum of two 
functions, one constant in the interval < 6 < 7r/M and the other piecewise constant 
and averaging to zero (these represent the sin and cosine fourier components defined 
in eq 5). The details will be given elsewhere; the results is that up to terms of order 
(ln(l/I)/I) 2 / 3 the energy is independent of the imposed phase O . 

This is related to the result of eq. (7), namely that for I 3> 1, the phase for 
any current of order Ij winds through may 2n revolutions between the center of the 
disk and the perimeter, so the average phase difference O is n °t an easily controlled 
quantity. To determine the critical current one must the energy subject to a constraint 
on the total current (which is no longer fixed by O ; one finds that the critical current 
is 2Ij/tt. 

Note that in the I > 1 case, unlike the I 1 solution with E < exists 

at 0o — 0. Therefore, as I is increased beyond some critical value, I c , presumably 
of order 1 the E < solution must appear. The transition region I ~ 1 requires 
numerical study, which has not been undertaken. However, by expanding eq. (5) to 
quadratic order is the c's and d's I have determined that the solution c n = d n = 
at 0o = becomes linearly unstable at I = I c = 1.53. Note also that at large I 
the 7r-periodicity which still exists in principle is of no practical relevance because if 
one starts the system in the lowest energy state at some initial phase O and then 
adiabatically increases the phase to 0o + 7r, the system will end up in a non-optimal 
state (higher than the optimum by a relative energy of order ln(I) 2 / 3 /I; to get to the 
optimal state would require overcoming an energy barrier of order Ija, which will be 
much larger than kt,T for a large device. 

I now extend the treatment of to include the effects of a finite London penetration 
depth, A, and of the energy stored in the magnetic field. One expects these effects to 
be important because currents in vortices decay exponentially on scales larger than 
the effective penetration depth A c fj (which for a film of thickness d < A may be greater 
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than the microscopic penetration depth A), and because the magnetic field extends 
also in the third dimension. The most interesting case is \ c s <C a; the currents in the 
vortices will then be confined to the region close to the A-B interface, so the precise 
geometry will not be important. It is convenient to study the linear interface shown 
in fig. 3. Here regions A and B are sheets of thickness d lying in the x-y plane and 
occupying the half planes x > and x < respectively. I assume the Josephson 
coupling Ij(y) along the A-B interface averages to zero, is periodic with period 2a, 
and satisfies Ij(y + a) = — Ij(y). To simplify the analysis I also assume that regions 
A and B have the same material parameters. 

I proceed as before, by first determining the general form of the bulk currents 
and then forcing this form to be consistent with the Josephson relation across the 
interface. In the presence of a magnetic field h one introduces a vector potential A 
given by y x A = h and one must use the gauge invariant supercurrent J = V0 — ^j^- 
instead of \j<p. It is also convenient to introduce a scaled field H via H = 27rh/$o 
Here <3>o — hc/2e is the superconducting flux quantum. In the volume occupied by 
the superconductor the equations governing the superflow and field are: V x J = H 
and v x H = J/A . Outside the film, J = and y x H = To simplify the analysis I 
assume that currents flow only in the x-y plane. This is equivalent to assuming that 
d > A or d < A. I give the analysis for the d > A case; the d <C A yields essentially 
identical results, but with A replaced by A e ff = A 2 /d. 

For d ^> A the solution in the film is only weakly perturbed by the spreading of the 

— * 

field in the region above and below the film. I therefore eliminate H from the in-film 

— * — * 

equations, solve for J, compute H in the plane of the film, and use this as a boundary 

— * 

condition for the equation for H in the region outside the film. The components of 

— * 

the current J are: 

n b p exp- Jl + b$\x\/\ 

J [x, y) = — > p , v 

r pity . v^y, , a v 

• [c p sin h dp cos J (8a) 



7j- i 

J v (x,y) = —sgn(x)^exp- yl + b 2 p \x\/\ 

(8b) 



. pny ixpy 
dp sin c v cos 



a a 
where b p = 7rpA/a. 

To implement the Josephson boundary condition one must compute the properly 
gauge invariant phase drop across the interface A0 GI . Now A0 GI has a contribution 
A x (Ax) from the change in the normal component of the vector potential across the 
interface. For a very thin junction, Ax — > and this contribution may be obtained 
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from the condition that the flux enclosed within the contour shown in fig. 3 is negligi- 
ble. (The usual contour extending to |x| = oo is not convenient because of the vortex 
currents). I find A0 GI (y 2 ) - A0 GI ( yi ) = J^ 2 dy J y (+Ax,y) - jr(-Ax,y). The free 
energy per 2a-period, normalized to the film thickness and to J c £o is given by the sum 
of three contributions: F = f curre nt + ffieid + fjos where f curre nt = / dxdy[J(x, y)] 2 , f fleid = 
>? j d 3 r H(f) 2 , and f Jos = I /* a ^g(y) cos[A(f) G i(6)} I have found that f fie i d < fcurrent, so 
I neglect f ne id henceforth. Using eqs. 8 and the Josephson boundary condition one 
sees that the free energy has the same form as that previously considered in eq. 5, 
except that the supercurrent term becomes 

instead of \ J2n n ( c n + d 2 )- 

Note that for b 2 ^> 1, eq. 9 reduces to the previous case, but for bp < 1, the 
supercurrent energy is much less than found previously. The previous analysis may 
now be applied with only minor modifications. The weak coupling limit previously 
meant I < 1; now the condition is I < A/a, which is much more stringent in high T c 
materials where A ~ 1400A and it is difficult to imagine deliberately made devices 
with scales a < 1/im ~ 10 4 A, although as discussed below wandering of the A- 
B interface may introduce a small length scale. In this weak coupling limit the 
maximum current is I 2 (a/ A), and the solution corresponds to vortices of extent A in 
the x direction and a in the y direction. In the range A/a < I < (a/A) a linearized 
equation cannot be used. Proceeding variationally one finds that modes with p < 
p max = (^pj have non-negligible amplitudes. (Note that b pmax < 1 if I < a/A). 
This implies that the scale of a vortex in the transverse (y) direction is ~ a/p max ~ 

M 1/2 „ f^ 1/2 



I / V It 



, i.e. is the Josephson penetration depth of the junction. The 
critical current is of order Ij Finally, for I > a/A the size of the vortex becomes less 
than the Josephson penetration depth and the problem becomes identical to the one 
previously considered; the solution corresponds to vortices of scale A in the transverse 
direction, pinched where they cross the A-B interface. 

I now turn to the experiment reported in ref. 1. For YBa 2 Cu 3 7 , £ — 10A [3] 
and Ij was measured [1] to be about 2 x 10 3 A/cm 2 , and the commonly accepted value 
is J c = 2 x 10 7 A/cm 2 . Because the films were apparently several thousand angrstoms 
thick, it seems reasonable to use the bulk value A ~ 2 x 10~ 5 cm for the penetration 
depth. The circumference, 27ra, of the device used in ref. 1, was 0.3 cm. Using these 
numbers I find I ~ 10 2 and a/A ~ 5 3 . It thus seems that the Josephson energy is 
dominant, but that the vortex size along the interface is controlled by the Josephson 
penetration depth Aj ~ 10~ 4 cm. Now in ref [1] "region A" was a hexagon and the 
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orientations were such that if the faces of the hexagon were flat on an atomic scale, 
then the Josephson coupling for two of the faces would vanish. The considerations 
presented here would then imply that would carry no supercurrent; the other four 
would carry a large current. In addition, the vortices at the corners of the hexagons 
should produce a magnetic field which extends a distance A away from the corner and 
Aj along the edge, and which at the surface would be of order $o/4AAj ~ lOgauss. 
In fact, the critical current through the individual faces was measured by destroying 
one face after another (via laser ablation) and monitoring the total critical current. 
The face-to-face variation of the critical current was found to be small, and for no 
face was the critical current found to be zero. This is not compatible with d-wave 
superconductivity if the interface is flat. However, if the interface wanders then the 
situation is different. In the rough interface case one might imagine that the Josephson 
coupling varies as one moves along the interface. If the length scale, L, over which 
the coupling varies is larger than the Josephson penetration depth then the previous 
arguments go though with the scale a replaced by L, and the experimental results 
might be consistent with d-wave superconductivity. 

Several extensions of the present work would be of use. A detailed examination 
of the hexagonal geometry used in ref. 1 should be performed. Also, a numerical 
solution in the region I ~ I c = 1.53 would be of interest. Finally, the dynamics 
of the system might be worth examining. If the system is forced to have a current 
greater than the critical current, voltage drops must occur. Presumably the bulk 
superconducting regions will be equipotentials, so there will be an angle-independent 
voltage drop across the A-B interface and the vortices shown in fig. 2 might for small 
voltages precess about the circle. 
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FIGURE CAPTIONS 



Fig. 1 System modelled in present paper shown for d-wave to d-wave case. Region A, 
inside the circle, is a film of d x 2_ y 2 superconductor with crystal axes running 
vertically and horizontally as shown Region B, outside the circle, is a film of 
d 2 . — y 2 superconductor with crystal axes rotated by 45°. The dots indicate 
places where the Josephson coupling across the circle vanishes by symmetry; 
the plus and minus symbols indicate the sign of the Josephson coupling in the 
regions between the dots. In the d-s case the Josephson coupling would vanish 
only at 0, ir and ±n/2. 

Fig. 2 Pattern of induced currents calculated at 0o — 7r /2 from eq. 9 for weak Joseph- 
son coupling for s-wave to d-wave case. Dots mark places where Josephson 
coupling vanishes. Inset shows motion of centers of vortices in presence of im- 
posed outward current. 

Fig. 3 Linear geometry used to study finite A case. Dots mark points at which a 
Josephson coupling changes sign. The +/— labels the local sign of Josephson 
coupling across interface. The dotted line is the contour used in the computation 
of A0GI- 
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